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We study several structure properties of finite nuclei using relativistic mean-field Lagrangians 
constructed according to the Brown-Rho scaling due to the chiral symmetry restoration at high 
densities. The models are consistent with current experimental constraints for the equations of 
' state of symmetric matter at both normal and supra-normal densities and of asymmetric matter at 

sub-saturation densities. It is shown that these models can successfully describe the binding energies 
, and charge radii of finite nuclei. Compared to calculations with usual relativistic mean-field models, 

these models give a reduced thickness of neutron skin in 208 Pb between 0.17 fm and 0.21 fm. The 
reduction of the predicted neutron skin thickness is found to be due to not only the softening of 
the symmetry energy but also the scaling property of p meson required by the partial restoration of 
chiral symmetry. 



o 
O 



(N 



PACS numbers: 21.10.-k, 21.60.Jz, 11.30.Rd 

I. INTRODUCTION 

The nuclear symmetry energy of isospin asymmetric nuclear matter not only plays a crucial role 
in a number of important issues in astrophysics, see, e.g., Rcfs. [1, 2], but is also important for 
■ understanding the structure of neutron- or proton-rich nuclei and the reaction dynamics of heavy- 
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ion collisions, see, e.g., Ref. [3-5]. However, the density dependence of the symmetry energy is still 



C ■ poorly known especially at supra- normal densities [6, 7]. In a previous work [8], we constructed 

several relativistic mean-field (RMF) Lagrangians using in-medium hadron properties according to 
f*^) ■ the Brown-Rho (BR) scaling due to the chiral symmetry restoration at high densities[9, 10]. The 



scalings and associated parameters that describe the in-medium hadron properties are consistent 
with those from microscopic calculations or those extracted from recent experimental data. The 
symmetric part of the resulting equations of state (EOS) around normal density is consistent with 
5^ ■ the data of nuclear giant monopole resonances [11] and at supra-normal densities it is constrained by 



the collective flow data from high energy heavy-ion reactions [12]. Moreover, the density dependence 
of the symmetry energy at sub-saturation densities is in agreement with that extracted from the 
recent isospin diffusion data from intermediate energy heavy-ion reactions [13-15]. An important 
feature of our models with the chiral limits is that the resulting symmetric nuclear matter EOS is 
soft at intermediate densities but stiff at high densities naturally, producing a maximum neutron 
star mass around 2.0Af Q consistent with the recent astrophysical observations. 

In the present work, we extend our well constrained RMF models to study ground-state properties 
of finite nuclei. In particular, we examine the neutron skin thickness of 208 Pb. The size of the 
neutron skin of finite nuclei is determined by the competition between the neutron pressure and 
the surface tension. Indeed, it has been shown in many studies that the neutron skin thickness 
of 208 Pb is rather sensitive to the density dependence of the symmetry energy [16-21]. The latter 
also influences the extraction of the incompressibility of nuclear matter from the experimental data 
of giant resonances in finite nuclei [19, 22, 23]. The size of the neutron-skin in 208 Pb predicted 
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by these models thus also appears to be limited by both the symmetry energy and the nuclear 
matter incompressibility of the models used. For instance, non-relativistic models that give an 
incompressibility of about k = 220 ± 15MeV and a symmetry energy at normal density less than 32 
MeV produce a neutron skin thickness for 208 Pb smaller than 0.2 fm, while the RMF models that 
have larger incompressibilities of 250-270 MeV and the symmetry energy at normal density larger 
than 34 MeV predict a neutron skin thickness for 208 Pb as large as about 0.3fm [21]. However, it is 
understood that it is the slope of the underlying symmetry energy that matters the most [21]. It was 
pointed out by Vretenar et.al. [19] that the incompressibility in RMF models was necessarily larger 
than 250 MeV for reproducing the nuclear structural properties. In a newly developed RMF model 
with the softened density dependence of the symmetry energy, Todd-Rutcl and Pickarcwicz obtained 
the incompressibility 230 MeV and predicted a neutron skin thickness of 0.21fm for 208 Pb [24]. On 
one hand, our models have the characteristic of a small symmetry energy of 31.6 MeV and a small 
incompressibility of 230McV at saturation density po = 0.16fm~ 3 . One would thus expect a thin 
neutron skin in 208 Pb. On the other hand, the cancellation of the medium effect from the vector 
meson mass and its coupling constant in our models results in an EOS in resemblance to those 
of the linear Walecka model and the best-fit model NL3 which have a much larger neutron skin 
thickness in 208 Pb. Given the successful description of the high-density behavior of nuclear matter, 
it is necessary to extend our models to the low density region by studying the ground-state properties 
of finite nuclei, especially the neutron skin thickness of heavy nuclei. Moreover, it is important to 
have a better understanding about the various contributing factors to the neutron skin thickness 
within the RMF models. 

The paper is organized as follows. In Section II, we briefly introduce the formalism of RMF with 
Brown-Rho scaling due to the chiral symmetry restoration at high densities. Results on the ground- 
state properties of finite nuclei, especially the neutron-skin thickness of finite nuclei are presented in 
Section III. A summary is finally given in Section IV. 

II. A BRIEF SUMMARY OF THE FORMALISM 

A lot of efforts had been devoted to describing well the binding energies and charge radii of finite 
nuclei simultaneously[25~28] within the RMF models. Inspired by the nonlinearity of a meson in 
chiral a model or induced by the model renormalization, nonlinear a meson terms were introduced 
to describe the medium effects that arc important in the Brucckncr theory. As a result, the binding 
energies and charge radii of finite nuclei were well reproduced, along with an improved description 
for the nuclear surface. Later on, the lu meson nonlinear term was introduced in compliance with 
the relativistic Brueckner results [29]. Some prominent representatives of the best- fit RMF models 
are NL3 [30], TM1 [29], FSUGold [24], and so on. Based on the microscopic relativistic Brucckncr 
results, the density-dependent Hartrcc or Hartrec-Fock approaches were developed to describe the 
finite nuclei (see [31, 32] and references therein). These approaches are usually parameter free but 
their predictions are not as accurate as the best-fit models. Our models used here have no nonlinear 
meson self-interactions, while the medium effects are given by the BR scaling. 

In the present work, the model Lagrangian with the density-dependent couplings and meson masses 
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is written as 

C = #7^ -M + ,g> - gl^uf - g* p l^ - e l -(l + r 3 ) 7 ^]V + \{d^a - mfa 2 ) 



7 V" + lm*Jco^ + \ m fb,X - - A A, U A^ (1) 



where ip,a,u, and bo are the fields of the nucleon, scalar, vector, and isovector-vector mesons, with 
their masses M, m*,m*, and m*, respectively. The meson coupling constants and masses with 
asterisks denote the density dependence, given by the BR scaling [8, 33, 34]. The Dirac equation of 
nucleons in the mean-field approximation is given by 

[i7 M 0" -M + g* a a - fl * 7o cJ - .g^o^o - e^(l + r 3 ) 7o ,4 + 7o £ fl <> = 0. (2) 

The rearrangement term E R, ° is essential for the thermodynamic consistency to derive the pressure 
and plays role in modifying the single particle energy and total binding energy of finite nuclei. In 
the mean field approximation, the is given by 



^ dC „ ,0ro* * h2 drn* p f 2 9m: dgZ dg* dg* a 



The density dependence of parameters is described by the scaling functions that are the ratios of 
the in-medium parameters to those in the free space. We take the scaling functions for the coupling 
constants of scalar and vector mesons as [8] 



1 * 1 -yp/po = ,s i-yp/p< 



1 + xp/po 1 + j/pp/ po 1 + VojP/Po 

For hadron masses, the scaling function reads 

Hp) = i - yp/po- (5) 

The energy per nucleon is given by [35] 

E = ^(2j a + l) Co - ^ f d 3 r(-g* a a Ps + g*u>p + g* p b p 3 + eA p p - 2^p) + E c . m . (6) 

where e is the eigen energy of nucleon, a denotes the occupied state, and j is the quantum number 
of total angular momentum. The energy from the center of mass correction is E c _ m , = —3/4 x 
41A -1 / 3 [28]. The equation of motion of mesons can be written out according to the Euler-Lagrange 
equation. The equations of nucleons and mesons are coupled nonlincarly and can be solved by 
iterations. The detailed procedure can be found in numerous literatures. Here, wc just emphasize 
that the meson mass in the Green function is that in the free space. As an example, the tr meson 
field is integrated out as follows 

a(r) = J dr x rl G(r,ri,m CT )[-g*(ri)/? s (ri) - (ml - ?< 2 (ri))cr(ri)] (7) 

where m a is the free a meson mass, g*(ri) = g*(p(ri)), and m* (ri) = m* (p(ri)). The expression of 
G can be found in Rcf . [361 ■ 



III. RESULTS AND DISCUSSIONS 



The details in determining model parameters can be found in our previous work [8] . Here we will 
firstly do calculations based on the parameter set SL1. In Table I, we tabulate the parameter set 
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TABLE I: Parameter sets used in the calculations. The vacuum hadron masses are M = 938MeV, 
m w = 783MeV and m p — 770MeV. The coupling constants given here are those at zero density. The 
incompressibility and the symmetry energy are respectively 230MeV and 31.6MeV at saturation density 
po = 0.16fm -3 for all models. The binding energy per nucleon (B/A), m a , and the incompressibility k are 
in unit of MeV. 



Model 


m a 


9° 




9p 


y 


X 


B/A M*/M 


y P y<*> 


SL1 


600 


10.3665 


10.4634 


3.7875 


0.126 


0.234 


-16.0 0.679 




SLC 


590 


10.1408 


10.3261 


3.8021 


0.126 


0.239 


-16.3 0.685 




SLCd 


590 


10.1408 


10.3261 


5.7758 


0.126 


0.239 


-16.3 0.685 


0.5191 - 



SL1 which is a little different from that in Ref. [8] because of a larger m a here. The nuclear matter 
property does not change at all by the variation of m a because the ratio ga/rria keeps unchanged, 
while the description for finite nuclei is modified. The latter can be seen in Eq. (7) and is correct 
for all cases irrespective of whether the coupling constant and meson mass are density dependent or 
not. For the linear Walecka model, the binding energies and charge radii can not be described well 
simultaneously by adjusting the m a [27]. The second term in the integrand of Eq. (7) is of prime 
importance in improving the description of the binding energies and charge radii of finite nuclei. 
Besides, another improvement comes from the inclusion of the rearrangement term appearing in 
Eqs. (2) and (6). Using m a = 500MeV in the original version of SL1, the total binding energy 
for 40 Ca is smaller than the experimental value by 70 MeV, and for 208 Pb the deviation is about 
230MeV. These results are improved appreciably by shifting the m a up to 600MeV, as shown in 
Table II. The relative error of charge radii is less than 1.5% and that of binding energy is less than 
3%. 

The accuracy of predictions can be much improved by setting B/A = — 16.3MeV and m„ = 
590MeV. Such a slight shift changes very little properties of infinite nuclear matter since the incom- 
pressibility is kept the same as that of the SL1 parameter set. The new parameter set called SLC 
is also listed in Table I. Except for 16 0, excellent agreement with the experimental data is achieved 
using the model SLC, as shown in Table II, though the binding energy of medium-mass nuclei is 
slightly underestimated by about 1.5%. To see the effects of the density dependence of the symmetry 
energy on the properties of finite nuclei, we derive the parameter set SLCd. The parameter set SLCd 
is different from the SLC only for the p meson coupling constant and its density dependence which 
is modified by an additional coefficient y p , as given in Table I. The results obtained by using the 
SLCd parameter set are given in Table II, where for comparison the FSUGold results are also listed. 
Fig. 1 shows the density dependence of the symmetry energy with the SLC, SLCd and SL1 param- 
eter sets. For comparison, we also include the MDI(x=0) and MDI(x=-l) results which represent 
the experimental constraints on the density dependence of the symmetry energy at sub-saturation 
densities, derived by studying the isospin diffusion data within a transport model [14, 15]. As shown 
in Fig. 1, no obvious difference between the SLC and SL1 is visible for the density dependence of 
the symmetry energy. Interestingly, all the symmetry energies from the parameter sets SLC, SLCd 
and SL1 are consistent with that constrained by the isospin diffusion data. 

It is now well known that the neutron skin thickness (r n — r p ) of heavy nuclei depends sensitively 
on the slope of the symmetry energy at normal density. A stiffer density dependence of the symmetry 
energy results in a thicker neutron skin. The difference between the symmetry energies with the SLC 
and SLCd, shown in Fig. 1, is responsible for the variation of the neutron skin in neutron-rich nuclei 
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FIG. 1: (Color online) The symmetry energy as a function of density for different models. The SL1 and SLC 
results overlap completely. 



TABLE II: The binding energy per nucleon (B/A), charge radii r c in fm, and neutron skins r n — r p obtained 
from different models. The available experimental data are from [37, 38]. 



Nucleus Expt. SL1 SLC SLCd FSUGold 



lb 


B/A 


7.98 


8.03 


8.07 


8.07 


7.96 




r c 


2.70 


2.72 


2.74 


2.74 


2.69 








-0.03 


-0.03 


-0.03 


-0.03 


40 Ca 


B/A 


8.55 


8.43 


8.54 


8.54 


8.54 




r c 


3.48 


3.43 


3.45 


3.44 


3.42 




T n Tp 




-0.05 


-0.05 


-0.05 


-0.05 


48 Ca 


B/A 


8.67 


8.41 


8.52 


8.46 


8.58 




r c 


3.47 


3.45 


3.46 


3.47 


3.44 




T n 'p 




0.20 


0.20 


0.18 


0.20 


90 Zr 


B/A 


8.71 


8.44 


8.59 


8.56 


8.68 




r c 


4.26 


4.22 


4.23 


4.23 


4.25 




T n 'p 




0.10 


0.10 


0.08 


0.09 


114 Sn 


B/A 


8.52 


8.24 


8.40 


8.38 


8.49 




r c 


4.61 


4.57 


4.58 


4.59 


4.60 




T*n Tp 




0.09 


0.09 


0.07 


0.09 


132 Sn 


B/A 


8.36 


8.18 


8.35 


8.23 


8.34 




r c 




4.68 


4.69 


4.70 


4.71 




i n Tp 




0.27 


0.27 


0.23 


0.27 


208 pb 


B/A 


7.87 


7.68 


7.87 


7.79 


7.89 




r c 


5.50 


5.46 


5.47 


5.48 


5.52 








0.21 


0.21 


0.17 


0.21 



as seen in Table II. Noting that the model FSUGold has almost the same symmetry energy in the 
whole density region as that of the MDI(x=0) and has larger symmetry energies at low densities than 
the SL1 and SLC, the neutron-skin thickness in 208 Pb with the SL1 and SLC would thus be expected 
to be much larger than that given by the FSUGold. However, from Table II, we can see surprisingly 
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that the SL1, SLC and FSUGold predict almost the same neutron skin for magic nuclei listed in 
Table II. Therefore, given the same density dependence of the symmetry energy as with the best-fit 
model FSUGold, the present models with the BR scaling predict smaller values for the neutron-skin 
thickness for finite nuclei. Compared to other best-fit RMF models, such as the NL1 and TM1 that 
predict the thickness of neutron skin in 208 Pb to be 0.28fm and 0.26fm, respectively, this reduction 
of the predicted neutron-skin thickness is more appreciable. To understand the underlying physics 
for the above observation, we may write down the integration form for the isovector potential that 
is a product of the p meson field and the coupling constant g* 

V P (r) =9*(r) J ^ 1 r?G(r,r 1 ,m p )[- 3 ;(r 1 )p 3 (r 1 )-(m2-m; 2 (r 1 ))6 (r 1 )]. (8) 

For the SL1 and SLC, the ratio C p = g* p /Tn* p is a constant. Properties of asymmetric nuclear matter 
in the RMF approximation depend on the constant ratio C p rather than the respective values of the 
coupling constant and the mass. This is the same as in the linear Walecka model and the best-fit 
model NL3. However, the linear Walecka model and the NL3 model predict much thicker neutron 
skin of 0.27fm and 0.28fm in 208 Pb, respectively [27, 30]. As the symmetry energy changes from 
35McV to 31.6McV in the linear Walecka model, the neutron thickness is still as large as 0.25fm. It 
is similar for the NL3 model, but a reduction in the symmetry energy to 31.6MeV results in large 
deviations of calculated nucleus masses from the experimental values [19], which is strongly unfavored 
by the best-fit model. The linear Walecka model, the NL3 and the FSUGold all have a much smaller 
effective nucleon mass M* at the saturation density than ours. It should be noted that it is not the 
larger M* that results in the small neutron skin since a small neutron skin of 0.22fm in 208 Pb is 
still obtained with our model that is adjusted to have a small M* = 0.6M with the same symmetry 
energies as given by the SL1 and SLC parameter sets. It is interesting that the scaling property of 
the p meson that is hidden behind calculations for nuclear matter is now displayed clearly in finite 
nuclei. The second term in the integrand of Eq. (8) is induced by the p mass scaling. Together 
with the density dependent coupling constant, it is responsible for the reduction of the neutron skin 
thickness in neutron-rich nuclei. The drop of the neutron skin thickness with the SLCd, compared to 
the SL1 and SLC, can just be attributed to the softening of the symmetry energy. We note that the 
neutron skin is totally insensitive to the value of m a , though a better fit of both binding energies and 
charge radii relies on the choice of m a . The thickness of neutron skin in 208 Pb varies from 0.21fm to 
0.17fm by changing from the SLC to the SLCd. This range is consistent with current measurements 
using the X-ray cascade from antiprotonic atoms: 0.16 ± 0.06fm [39] and also agrees well with that 
obtained from the analysis of the isospin diffusion data [20]: 0.22 ± 0.04. Interestingly, the predicted 
range of neutron skin in 208 Pb between 0.21fm and 0.17fm is very close to that predicted by using 
some Skyrme interactions [21]. 

It is worth noting that a much smaller neutron-skin thickness in 208 Pb can be obtained within 
the RMF models provided that smaller values of symmetry energy are given. For instance, the 
ES25 model [2] constructed with the incompressibility of 211.7 MeV and the symmetry energy of 
25 MeV at saturation density predicts a neutron skin as thin as 0.138 fm for 208 Pb. With the same 
incompressibility and symmetry energies as the ES25 model, our model predicts an even smaller 
neutron skin of 0.126 fm for 208 Pb. Of course, the magnitude of the reduction due to the in-medium 
properties of the p meson is accordingly smaller. 

The deficiency of our models mainly lies in the unsatisfactory description for spin-orbit splittings. 
For instance, the lp spin-orbit splitting in 16 is 3.7 and 3.5MeV with the SL1 and SLC, respectively, 
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which is smaller than the experimental value by about 2.5 MeV. This deficiency stems mainly from 
a much larger Dirac mass M* in our present models, i.e., about 0.68M. Usually, a Dirac mass of 
about 0.5-0.6M can reasonably describe the spin-orbit splitting in RMF models (see, e.g., Ref. [40]). 
As a result, the spin-orbit potential of the central part is much suppressed in our models, compared 
to that given by the best-fit models. The deviation of spin-orbit splittings from the available data is 
comparable to that given by the density dependent approaches based on the relativistic Brueckner 
results [32] . The right spin-orbit splitting in our calculations can be obtained by reducing the Dirac 
mass moderately. On the other hand, the larger Dirac mass leads to larger Landau mass (about 0.74 
in our models), which gives a more reasonable description for the level density [40]. 

IV. SUMMARY 

In summary, we have studied ground-state properties of finite nuclei within the relativistic mean- 
field models that are constructed according to the Brown-Rho scaling due to the chiral symmetry 
restoration at high densities. Not only can the constructed models SLC and SLCd be consistent 
with current experimental results for the equation of state of symmetric matter at normal and 
supra-normal densities and of asymmetric matter at sub-saturation densities, but also present a 
fairly satisfactory description for the ground state properties of finite nuclei. The binding energies 
and charge radii of a variety of magic nuclei are nicely reproduced. Our results indicate that the 
neutron-skin thickness of finite nuclei depends not only on the density dependence of the symmetry 
energy but also on the scaling property of the p meson. Compared to calculations of usual RMF 
models a reduction of neutron skin thickness of neutron-rich nuclei is observed in our models. We 
find that the scaling property of the p meson that is hidden in nuclear matter calculations is shown 
up clearly by the reduction of the neutron skin in 208 Pb. The neutron skin thickness of 208 Pb is 
predicted to be between 0.17 fm to 0.21 fm consistent with the current measurements. 
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